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ABSTRACT 



The observational study of star formation relations in galaxies is central to unraveling the related physical processes that are at work 
on both local and global scales. It is still debated whether star formation can be described by a universal law that remains valid 
in different populations of galaxies. We wish to expand the sample of extreme starbursts, represented by local luminous and ultra- 

O luminous infrared galaxies (LIRGs and ULIRGs), with high quality observations in the I~0 line of HCN, which is taken as a proxy for 
the dense molecular gas content. The new data presented in this work allow us to enlarge in particular the number of LIRGs studied in 
, HCN by a factor 3 compared to previous works. The chosen LIRG sample has a range of HCN luminosities that partly overlaps with 

that of the normal galaxy population. We study if a universal law can account for the star formation relations observed for the dense 
molecular gas in normal star forming galaxies and extreme starbursts and explore the validity of different theoretical prescriptions 
of the star formation law. We have used the IRAM 30m telescope to observe a sample of 19 LIRGs in the 1-0 lines of CO, HCN 
^ , and HCO*. The galaxies have been extracted from a sample of local LIRGs with available high-quality and high-resolution images 

5_( ■ obtained at optical, near and mid IR wavelengths, which probe the star formation activity. We have thus derived the star formation 

' rates using different tracers and determined the sizes of the star forming regions in all the targets. 

' The analysis of the new data proves that the efficiency of star formation in the dense molecular gas (SFEdcnsc) of extreme starbursts 

is a factor 3-4 higher compared to normal galaxies. Kennicutt-Schmidt (KS) power laws have also been derived. We find a duality in 
KS laws that is further reinforced if we account for the likely different conversion factor for HCN (a'^''^) in extreme starbursts and 
for the unobscured star formation rate in normal galaxies. This result extends to the higher molecular densities probed by HCN lines 
^ ' the more extreme bimodal behavior of star formation laws, derived from CO molecular lines by two recent surveys. 

CO , We have confronted our observations with the predictions of theoretical models in which the efficiency of star formation is determined 

. by the ratio of a constant star formation rate per free-fall time (SFRff) to the local free-fall time (tfj). We find that it is possible to fit 

' the observed differences in the SFEjensc between normal galaxies and LIRGs/ULIRGs using a common constant SFRjf and a set of 

' physically acceptable HCN densities, but only if SFRff~0.005-0.0I and/or if a^™ is a factor of ~a few lower than our favored values. 

Star formation recipes that explicitly depend on the galaxy global dynamical time scales do not significantly improve the fit to the 
new HCN data presented in this work. 

Key words. Galaxies: evolution - Galaxies: ISM - Galaxies: starburst - Infrared: galaxies - ISM: molecules - Radio lines: galaxies 
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1 . Introduction (Pgas): Psfr P^as- Schmidt ( 1 19591 1 predicted that should be 

around 2. If we suppose that the relevant time scale for star 
The quest for a quantitative description of how the star formation formation is the local free-fall time (tff), in this case psFR 
rate should depend on a set of galaxy parameters (gas and stellar ^ Qn the condition that the gas scale-height is con- 
content, gas kmematics, environment) is a fundamental though ^^^^^ ^^^^^ galaxies, a non-trivial additional assumption, the 
not fully understood problem. Observational constraints ai-e cen- equivalent of the Schmidt law in terms of the corresponding sur- 
tral to understanding the role that small-scale (local) and large- f^^g densities of SFR (Ssfr) and gas (S„as) is Ssfr <^ S'f,. 
scale (global) galaxy properties play in star formation relations. ^ 
This is a requisite to validating the recipes of star formation used Krumholz &McKee C2005) (see also Krumholz, McKee & 
in numerical simulations of galaxy evolution. Tumlinson 2009 ) developed a model that successfully describes 
Based on simple theoretical grounds Schmidt (119591 1 pro- how star formation takes place in a highly turbulent molecular 
posed that the star formation rate (SFR) per unit volume (psfr) medium. Their scenario pictures star formation as a very ineffi- 
should be a power law of index A^ of the gas volume density cient process that takes place in small subregions of supersonic 

turbulent and virialized molecular clouds. The surface density 

* Based on observations caiTied out with the IRAM 30m telescope, of the star formation rate is in any case expected to scale as a 

IRAM is supported by INSU/CNRS (France), MPG (Germany) and power law of the mean density of the gas with an index close 

IGN (Spain) to 1.5. Star formation in Krumholz &McKee's formulation is 
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regulated by processes mostly within the clouds. The relevant 
time scale for these models, the free-fall time for cloud col- 
lapse, is thus a local time scale. In this type of formulation, 
star formation laws should not depend explicitly on global dy- 
namical time scales or be determined by galactic-scale phenom- 
ena. Instead, large-scale processes are thought to be responsible 
for making the gas assembly prior to star formation in galax- 
ies. Alternatively, other forms of the star formation laws ad- 
duce the importance of including expUcitly global dynamical 
time scales tdyn ( Elmegreen [19971 Silk (19971 Kennicutt [19981 
Elmegreen llofel Boiss ier et al. [20031 Tasker & Tan [20091 
Tan 12010 1 Genzel et al. I2010I I. The star formation rate density 
SsFR should then scale as o^Sgas/tdyn- It is yet unclear if the two 
forms of the star formation laws are fully equivalent both in 
normal galaxies and in mergers. Kennicutt (1998 ) found that if 
tdyn is the average orbital dynamical time scale, both formula- 
tions are equally successful at describing star formation laws in 
the disks of normal galaxies. The physical framework for this 
equivalence relies on the fact that normal galaxy disks, viewed 
as dynamically relaxed stable systems, fulfill the Toomre sta- 
bility criterium (Toomre [l964| l. It is nevertheless questionable if 
this equivalence should also hold in extreme starbursts which are 
often hosted by merging systems. 

During the last decades observers have been using differ- 
ent proxies for Ssfr and Egas to validate the predictions of 
star formation law models. These observations have been con- 
ducted on different galaxy samples, including normal galaxies 
and extreme starbursts at different redshifts. They also encom- 
pass a wide range of spatial resolutions, which allow not only 
global measurements but also studies of spatially resolved ob- 
jects. The different observers do not find a single exponent fit- 
ting the entire star formation relation. Kennicutt d 19981 1 finds 
an index n= 1.40+0. 15 using CO and HI as gas tracers and a 
~100 galaxy sample that includes normal galaxies, nuclear star- 
bursts and a few LIRGs (~5) at z=0. Yao et al. ( 12003]) find an 
index n=1.40±0.30 in their CO survey of IR luminous galax- 
ies. Bouche et al. (120071) . by including submillimeter galax- 
ies (SMGs) and using a lower CO conversion factor in z=0 
ULIRGs, determines a higher index for the Kennicutt-Schmidt 
(KS) law n~1.7. The high-resolution survey of 18 normal galax- 
ies published by Bigiel et al. (lOOSi and Leroy et al. (I2008I I 
found a linear KS relation above a gas surface density thresh- 
old of IOMqpc"^, which determines the transition from atomic 
to molecular gas, here traced by CO lines. 

More recently the CO surveys of Daddi et al. ( 120101 ) and 
Genzel et al. (1201 Oi l presented evidence that normal galaxies 
and mergers (LIRGs, ULIRGs and SMGs) occupy different re- 
gions in the molecular gas mass versus star formation rate plane. 
These results suggest the existence of a bimodality in star for- 
mation laws, where normal galaxies show a 4-10 longer deple- 
tion time scales compared to mergers. While the KS relation 
remains mostly linear inside each galaxy population, bimodal- 
ity introduces a discontinuity in the two-function power law. 
These authors discussed that when global dynamical time scales 
are included a universal star formation law is obtained. Based 
on a study of far-infrared fine structure lines, Gracia-Carpio et 
al. |201 1) concluded that extreme starbursts display significantly 
enhanced line deficits. The results of Gracia-Carpio et al. (201 Jj 
provide additional support for the existence of different proper- 
ties of the interstellar medium in normal star forming galaxies 
and extreme starburst systems. 

Krumholz & Thompson (|20071 l (see also Narayanan et 
al. [2008 ) studied how the power index of KS laws determined 
from observations should change depending on critical density 



of the tracer used to probe the star forming gas, a prediction 
confirmed by observations done in several dense gas tracers 
(Narayanan et al. [20051 Bussmann et al. 2008 ; Bayet et al. 20091 
Juneau et al. 120091 ). In the particular case of HCN(l-O), only in 
galaxies where the average gas density exceeds a few 10"* cm"^ 
(i.e., the effective critical density of the HCN J=l-0 line), we 
would start recovering the expected superlinear behavior of the 
universal KS law derived by Krumholz & McKee (2005 ). 

Observations of HCN lines have been used to study dense 
molecular gas properties in samples of normal spirals, LIRGs, 
ULIRGs, QSOs and molecular gas-rich high-z galaxies (e.g., 
Gao & Solomon [2004ai2004bl Solomon & Vanden Bout [20^ 
Wagg et al. 200T, Riechers et al. [2007l Gracia-Carpio et al. [20061 
2008). By contrast to CO line observations, which are sensitive 
to the global molecular gas content, HCN observations are a 
much better tracer of the dense molecular phase, which is more 
directly related to star formation, and can thus put stringent con- 
straints on star formation models. The most important result pre- 
sented by Gao & Solomon in their seminal papers is the discov- 
ery of a tight correlation between the infrared and the HCN(1- 
0) luminosities over three orders of magnitude in Lir. This result 
was interpreted as evidence that star formation is the main power 
source in ULIRGs. The linearity of the tight coiTelation between 
the IR and the HCN(l-O) luminosities implies that the star for- 
mation efficiency measured with respect to the dense molecular 
gas content (SFEdense <^ ^ir/^hcn(i-0)) constant in all galax- 
ies, independently of Lir. 

More recent results have partly questioned this picture, how- 
ever. Gracia-Carpio et al. (2008) published observations made 
with the IRAM 30m telescope of the J=l-0 and 3-2 lines of 
HCN and HCO^ used to probe the dense molecular gas content 
of a sample of 17 LIRGs and ULIRGs. These observations were 
also used to derive a new version of the power law describing 
the correlation between Lfir and i[jcN(i-()) f™!^ normal galax- 
ies (Lir < 10' ' Lq) to high-z galaxies. The results of Gracia- 
Carpio et al. (2008) indicated that the Lfir/L^j^j^^-j ratio, taken 
as proxy for SFEdense, is a factor ~2-3 higher in IR luminous tar- 
gets compared to normal galaxies. Furthermore, Gracia-Carpio 
et al. ( 120081) found, based on a multiline analysis of HCN and 
HCO^ data, that Xhcn is probably about 3 times lower at high 
Lfir- Taken together, these findings suggest that SFEdense may be 
up to an order of magnitude higher in extreme starbursts (LIRGs, 
ULIRGs) than in normal galaxies. 

Compared to normal galaxies and the most extreme star- 
bursts represented by ULIRGs, LIRGs can constitute the tran- 
sition point in the star formation laws of galaxies. The main pur- 
pose of this paper is to enlarge the sample of LIRGs with high 
quality HCN(l-O) observations. This is a key to testing if the 
observed bimodality of star formation laws in galaxies derived 
from CO, can be extended to the higher density regime probed 
by HCN. We use these results to explore the validity of different 
star formation law theoretical prescriptions which are cuiTently 
debated. To this aim we have used the 30m telescope to observe 
a sample of 19 LIRGs in the 1-0 fines of CO, HCN and HCO+. 
Most of the targets have been extracted from the complete sam- 
ple of local (d<78 Mpc) LIRGs studied by Alonso-Herrero et 
al. (!2006). All the targets in this sample have high-quality and 
high-resolution imaging at optical, NIR and MIR wavelengths 
(obtained with HST, VLT, CAHA and Spitzer). This allows us 
to derive accurately the sizes of the star forming regions and the 
star formation rates for all the targets. The new data presented 
in this work allow us to expand the number of LIRGs studied in 
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Table 1. Main properties of the new sample of LIRGs. 
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Notes. " t/sF is the source size (major axis) of the star forming region in arscec determined from high-resolution optical/NIR images of hydrogen 
recombination lines (Ha, Paa). 



HCN-HCO^ lines by more than a factor 3 compared to previous 
works (9 LIRGs in Gracia-Carpio et al.'s sample). 

2. The sample of LIRGs 

We have obtained new observations for a sample of 19 LIRGs. 
Out of these, fourteen have been extracted from the sample of 
nearby LIRGs studied by Alonso-Herrero et al. (120061 ). This is 
a volume-limited sample (v<5200 km/s, distances of between 
35 and 78 Mpc) drawn from the IRAS Revised Bright Galaxy 
Sample (RBGS) analyzed by Sanders et al. (2003). Alonso- 
Herrero et al. (2006) required the logarithm of the total IR lu- 
minosity to be >11.05 L0 and the galaxies to be at Galactic lati- 
tude b> 10°. The remaining five LIRGs were selected to populate 
the high end of the IR luminosity distribution of LIRGs. For the 
Alonso-Herrero et al. LIRGs we have obtain ed HST/NICMOS 
Paa imaging (Alonso-HeiTero et al. 120011 '2002, 2006), as 
well as ground-based optical integral field spectroscopy using 
CAHA/PMAS and VLT/VIMOS covering the Ha spec tral re - 
gion (Alonso-Herrero et al. 2009; Rodriguez Zaurm et al. l201 II ). 
In addition there are mid infrared (MIR) Spitzer imaging obser- 
vations for the entire sample from the GOALS sample (Armus et 
al. 120091 ) ■ This means that the star formation rates (SFR) of the 
sample are very well characterized from hydrogen recombina- 
tion lines (Ha, Paa) and MIR (24jum) luminosities on the same 
angular scales probed by the IRAM observations. We also have 
a good estimate of the typical size of the star forming regions 
in all the galaxies of our sample, thanks to the available high- 
resolution imaging (Ha, Paa). This minimizes the errors when 
both SsFR and Sdense, required to derive KS laws, are evaluated. 
The estimated sizes are typically <10" in all sources, equivalent 
to 1.7-3.6 kpc for the given range of distances. This guarantees 
that the total emission of the molecular gas is contained in the 
30m beam (-28" at the frequency of the HCN(l-O) line). 

Additionally, in this paper we use a compiled sample of 
108 galaxies with published FIR and HCN(l-O) observations 
in the literature. When available, we also include HCO^(l-O) 



data. Besides the new data obtained in LIRGs for this work 
we include the data of normal galaxies, LIRGs and ULIRGs 
of Gao & Solomon (I2004al ). data from the sample of infrared- 
excess Palomar-Green QSOs (Evans et al. 2006) and data from 
high-z galaxies with available HCN observations (see Gao 
et al 120071 and references therein). Following the same ap- 
proach adopted by Gracia-Carpio et al. (I2008I I. in those high-z 
sources where HCN(l-O) are not available, ^hcnci-O) 
rived assuming that their rotational line luminosity ratios are 
^HCN(2-i)/^HCN (i-0) ~ similar to the mean value measured 
by Krips et al. mm . For the z ^ 4 quasar APM 08279+5255 
we adopt ^hcn(5-4)/''"hcn(1-0) ~ assuming the physical con- 
ditions derived in Garcfa-Burillo et al. (12006b . 

3. Observations 

The new molecular line observations of the sample of 19 LIRGs 
described in Sect. |2] were carried out in one observing run on 
June 2008 with the IRAM 30m telescope at Pico de Veleta 
(Spain). We tuned the 3 mm and 1 mm SIS receivers of the 
30m telescope to the redshifted frequencies of the HCN(l-O), 
HCO+(1-0), CO(l-O) and CO(2-l) lines. The beam sizes of the 
30m telescope range from FWHM~28" at 88 GHz to -11" at 
230 GHz). We covered a velocity range 1300 - 1800kms ' for 
the 3 mm lines and 1200km s"' for the 1 mm lines. We used the 
wobbler switching mode to make baselines flatter System tem- 
peratures during the observations were typically ~100-130K at 
3 mm and ~300-500K at 1 mm. Receivers were used in single 
side-band mode (SSB), with a high rejection of the image band: 
>12dB at 1 mm and >20dB at 3 mm. The calibration accuracy 
is estimated to be better than 20%. We checked the pointing of 
the telescope every 1.5 hours by observing nearby continuum 
sources; the average rms pointing error was 2" during the ob- 
serving run. 

Throughout the paper, velocity-integrated line intensities (/) 
are given in antenna temperature scale, T*. The T* scale re- 
lates to the main beam temperature scale, Tmb, by the equa- 
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Fig. 1. a) CO(l-O) spectra observed with the IRAM 30m telescope in a sample of 19 LIRGs. Spectra appear from top to bottom 
and from left to right in order of increasing right ascension of the source. To ease the display some of the line intensities have been 
scaled by the factor indicated in the panel. Grey-filled histograms highlight the velocity range used to calculate the baseline fitting 
and the line areas. 



tion Tmb = (Fe,ff/Bfis)T*, where F^ff and B^s are the forward 
and beam efficiencies of the telescope at a given frequency. For 
the IRAM 30m telescope Feif/Beif = 1-22 (1.57) at 87-88 GHz 
(230GHz) and S/T^^ = 4.95 JyK The velocity windows 
used to derive /hco+ and /hcn have been defined from the corre- 
sponding higher signal-to-noise CO line profiles. We have sub- 
tracted linear baselines to the individual scans and eliminated 
those showing instabilities requiring higher order baselines. We 
show in Fig. [T]the ensemble of HCN(l-O), HCO+(1-0) and 
CO(l-O) spectra used in this work. Molecular line luminosities 
(L') were computed in units of L' - Kkms 'pc^ as defined 
by Gao & Solomon (I2004al l. Luminosity distances have been 
derived assuming a flat A-dominated cosmology described by 
Ho = 71kms ' Mpc ' and Q.^ = 0.27. Observational parame- 
ters and results are summarized in Tables[T]and|2] 



4. Star formation laws of the dense gas 



In Sects 14.11 and 14.21 we study the star formation efliciency 
of the dense gas (SFEdense) and the Kennicutt-Schmidt (KS) 



laws in the different populations of galaxies compiled for this 
work, including the new sample of LIRGs. Following the 
same approach adopted by Gracia-Carpio et al. ( I2008I I. we de- 
rive SFR [M0yr-i] from L^^ defined as LfP'*=1.3x Lfir(40- 
500 yum). This aims at minimizing the possible AGN contribu- 
tion to LiR(8-1000/im). The AGN contribution to Lir in the 
population of local LIRGs with AGN detections is nonetheless 
known to be small, with a median value of Lboi(AGN)/LiR=0.05 
(Alonso-Herrero et al. |201 IK 

In this section we adopt the common assumptions regard- 
ing conversion factors used to derive SFR from and Mjense 
from ^HCN (i-O)- ^^^^ adopt the same universal factor used by 
Kennicutt dT998] l: 



SFR[Moyr-'] 



1.7xlO-«'LfR^[Lo] 



(1) 



Similarly, Mdense is obtained in the following from incNd-O) 
assuming a universal conversion factor (o-*^'--'^) similar to that of 
Gao & Solomon (.2004aJ : 



ense 

[Mq] - 10L[jQ,j(i_()) [Kkms ' pc ] 



(2) 
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Fig. 1. b) Same as a) but showing the HCN(l-O) spectra. 



We nevertheless question the grounds for these two working 
assumptions and the consequences thereof in Sect. |6] 



■^HCN(l-O) 



^ 0.59 D 



0.29 



(4) 



4.1. Star formation efficiency fSFEdensej 

Figure 12] shows the Lpij^/ L'^^^j^^^ ^^ luminosity ratio versus Lfir 
for the different populations of galaxies in this compilation, de- 
fined in Sect. |2] We adopt here the common assumption that 
the LpiR/iHCNCi-O) luminosity ratio is a good proxy for the star 
formation efficiency of the dense gas (SFEdense=SFR/Mciense 
^'Hr/^'hcnci-O))- estimate that the total uncertainty on both 
Lfir and ^^^^^(1-0) ^'^^'^ related quantities), with an equal 
contribution of statistical and systematic errors, amounts to 
±30% (±0.13 dex). The corresponding total uncertainty on 
Lir/l;^cn(i-0) is ±42% (±0.19 dex). 

In Fig.|2]we can see that SFEdense increases with Lfir from 
normal galaxies to extreme starbursts, represented by LIRGs, 
ULIRGs and high-z objects. In particular, the orthogonal regres- 
sion fit to the full sample of objects gives a power law index for 
SFEdense significantly different to zero: 



log 



-^FIR 



V HCN(l-O) J 



(0.29 ± 0.04) log Lfir - (0.23 + 0.42) (3) 



The location of the new sample of LIRGs in the SFEdense di- 
agrams of Figs. |2]and[3j' coiToborates on a more solid statistical 
basis, compared to the results of Gracia-Carpio et al. (2008), that 
SFEdense IS on average a factor ~2-3 higher in LIRGs/ULIRGs 
compared to normal galaxies (<SFEdense>(LIRGs/ULIRGs) 
~1400±100LoL'"'; <SFEdense>(normal)~600±70 LqL'"'). 
This difference is a factor of 2 larger than the total uncertainty 
on individual data points in Fig. |2] The total change in 
SFEdense cxpands an order of magnitude from normal galaxies 
to LIRGs/ULIRGs. Of particular note, two objects of the new 
sample, NGC 1614 and II Zw96, show the highest SFEdense val- 
ues (>3500LoL' ') thus far reported in a galaxy. These values 
are close the upper limit imposed by the maximum efficiencies 
measured in Galactic star-forming cores (~3900LoL'"'; Wu 
et al. I2005I I. The resulting index derived for the SFEdense law 
(0.29±0.04) is larger than the one obtained by Gracia-Carpio et 
al. (120081) (0.24±0.04). 
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Fig. 1. c) Same as a) but showing the HCO^(l-O) spectra. 



Figure [3^ shows the LpiR-iecNii-O) correlation plot obtained 
for the full sample of galaxies. The orthogonal fit to the data 
points gives a super-Unear correlation: 



logLpiR = (1.23 + 0.05) log L^jcNd-o) + (1-07 + 0.40) 



or L] 



FIR 



12 L' 



1.23 

HCN(l-O) 



(5) 



(6) 



The fit is similar to that found by Gracia-Carpio et al. (I2008I I. 
With the addition of the new LIRG sample we thus confirm that 
the LFiR-i'HCN(i-O) correlation is significantly super-linear. This 
partly contradicts the first findings of Gao & Solomon ( i2004al 
I2004bl l who determined a power index about 1 for the Lfir" 
^HCN(i-O) scatter plot and a correspondingly constant SFEdense 
independent of galaxy type. 

Figure |4]shows the LpiR-i^Hco+ci-o) correlation plot obtained 
for the sample of 34 LIRGs/ULIRGs where we have obtained 
data in the HCO^ lin43 (data from this work and from Gracia- 



' There is no equivalent survey of HCO^ data in a significant number 
of normal galaxies. The comparison between the star formation proper- 
ties of the dense molecular in LIRGs/ULIRGs and normal galaxies has 
thus to rely hereafter mainly on HCN data. 



Carpio et al. '2006', '2008). The orthogonal fit to the data points 
gives a close to linear correlation: 

logLpiR = (1.06 ± 0.15) log l;jco-(i_o) + (2.64 + 0.04) (7) 



or L] 



FIR 



: 436 L' 



1.06 

HCO+(1-0) 



(8) 



We note that a similar fit using HCN(l-O) data but restricted 
to the sample of LIRGs/ULIRGs gives also a close to linear 
power law (n=1.10±0.10). A similar fit of the Lfir-^^hcnci-O) 
relation restricted to normal galaxies also gives a power index 
compatible with unity: «=0.95±0.20. Altogether this is an in- 
dication that it is the inclusion of both normal galaxies and 
LIRGs/ULIRGs in the fit of Eq. |5] that makes the relation be- 
tween LpiR and i'HCN(i-o) become superlinear. 

Figure |5] shows the dense molecular gas depletion 
or exhaustion time scale tdep in Myr units, defined as 
tdep=SFEjJ^^^=Mdense [MqI/SFR [MQyr-'], as a function of 
LpiR. NGC 1614 and II Zw96 lie at the lower end of the tdep dis- 
tribution and are characterized by significantly short depletion 
time scales ^10 Myrs. This is an indication that they represent 
extreme starburst systems that will exhaust their dense molecular 
gas content on time scales comparable to the typical dynamical 
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Table 2. Results derived from the IRAM 30m telescope observations. 



galaxy 


Ico 


a 


'HCN 


Jhco* 




Kkms-' 


Kkms-' 


Kkms-' 


NGC 23 


28.36 


(0.32) 


0.81 


(0.13) 


1.08 


(0.13) 


MCG+ 12-02-001 


36.50 


(0.39) 


0.89 


(0.10) 


1.34 


(0.10) 


III Zw 35 


10.43 


(0.35) 


1.08 


(0.15) 


0.50 


(0.10) 


UGC 1845 


32.90 


(0.31) 


1.29 


(0.13) 


1.00 


(0.13) 


NGC 1614 


27.16 


(0.47) 


0.79 


(0.11) 


1.56 


(0.11) 


UGC 3351 


46.51 


(0.40) 


2.17 


(0.13) 


1.81 


(0.13) 


NGC 2388 


42.60 


(0.52) 


1.51 


(0.11) 


1.69 


(0.11) 


MCG+02-20-003 


14.08 


(0.30) 


0.52 


(0.10) 


0.54 


(0.10) 


NGC 2623 


19.02 


(0 34) 


1.30 


(0 11) 


0.67 


fO 12) 


NGC 3110 


31.37 


(0.46) 


0.92 


(0.10) 


0.87 


(0.10) 


IRAS 10173+0828 


3.52 


(0.26) 


0.43 


(0.11) 


0.29 


(0.11) 


TP CAA 
iV^ ouu 


9.61 


(0.25) 


0.80 


(0.10) 


0.55 


(0.10) 


NGC 5653 


36.56 


(0.57) 


1.00 


(0.10) 


1.00 


(0.10) 


NGC 5936 


27.44 


(0.30) 


0.94 


(0.08) 


0.73 


(0.10) 


IRAS 17138-1017 


27.36 


(0.35) 


0.58 


(0.09) 


1.12 


(0.13) 


IRAS 18090+0130 


19.06 


(0.34) 


0.51 


(0.11) 






NGC 6701 


36.17 


(0.32) 


1.11 


(0.17) 


1.15 


(0.08) 


II Zw 96 


8.97 


(0.32) 


0.35 


(0.09) 






NGC 7591 


24.66 


(0.31) 


0.61 


(0.13) 


0.57 


(0.13) 



Notes. " All velocity integrated intensities are in T* scale. We list l-cr errors within brackets. 




Lfir [Lq] 



Fig. 2. a) (left panel) We plot the Lfir/L^^j^^j luminosity ratio, a proxy for the SFEjense defined in Sect.|2l as a function of LpiR 
for different galaxy samples compiled in this work (See Sect. |2] for details). The different symbols represent the normal galaxies 
and L IRGs/ ULIRGs published by Gao & Solomon (2004a 2004b) (crosses), the LIRGs/ULIRGs published by Gracia-Carpio 
et al. (I2OO8I 1 (empty circles), PG QSOs and high-z galaxies published and compiled by Evans et al. (12006b and Gao et al. (2007ll, 
respectively (triangles), and the new sample of LIRGs published in this work (filled circles). Arrows represent upper and lower limits 
to ^HCN(i-o)- {fight panel) Same as a), but limits have not been represented. The solid line visualizes the orthogonal regression 
fit calculated for the full sample of objects, not taking into account limits. We highlight the location of NGC 1614 and II Zw96 in 
this plot. The total uncertainty of individual data points, including statistical and estimated systematic errors, are indicated by error 
bars: ±0.13 dex(±30%) in LpiR and ±0.19 dex(±42%) in Lpir/L'^-oy 



time scales for these sources (see discussion in Sect. 17.31 and in 
Sect. 3.6 of Alonso-Herrero et al.ilOOljfor NGC 1614). 



4.2. Kennicutt-Schmidt power laws 

Figure|6]shows the KS law obtained for the different populations 
of galaxies compiled for this work., i.e., the star formation rate 
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Fig. 3. a) (left panel) LFiR-Lyf,j^([_Qj correlation plot derived for different galaxy samples. The solid line visualizes the orthogonal 
regression fit to the full sample showing a super-linear correlation. The errorbars show the total uncertainty of individual data points 
(+0.13 dex=+30% in both Lhr and L'^qj^) ■ b) (right panel) We plot the Lpiji/ L'^^^j^i^^ ^^ luminosity ratio as a function of i^jov,. 
We display with different symbols normal galaxies (Lir < 1O"L0, open squares) and luminous infrared galaxies (Lir > 1O"L0, 
filled squares). The location of the new sample of LIRGs in this diagram is identified by blue color markers. The dashed and 
continuous horizontal lines indicate, respectively, the average value of Lfir/L^^j^^j in normal galaxies (~6OO+7OL0L' ') and 
LIRGs/ ULIRGs (-UOO+IOOLqL' '). The errorbars show the total uncertainty of individual data points: ±0.13 dex(±30%) in 
L^cN ±0.19 dex(±42%) in LFiR/^HCNd-O)- 
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Fig. 4. LpiR-i'Hco+d-o) coiTelation plot derived for a subset of 
LIRGs/ULIRGs. The solid line visualizes the orthogonal regres- 
sion fit to the sample. Errorbars as in Fig|3] 
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Fig. 5. The dense molecular gas depletion or exhaustion time 
scale, tdep=l/SFEdense in Myrs, derived for normal galaxies and 
LIRGs/ULIRGs. Symbols are as in Fig.|3l We highlight the loca- 
tion of NGC 1614 and II Zw 96 in the diagram. Errorbars show 
total uncertainties: +42% in tjep and +30% in Lfir. 



surface density Ssfr as a function of dense molecular gas surface 

density as traced by the HCN(l-O) line, Xdense- In order to derive published CO or HCN interferometer maps, which are available 
SsFR and Sdense wc adopted the molecular gas size estimates from for most of the sources culled from the literature. As for the new 
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sample of LIRGs, we estimated the size of the central star form- 
ing region from the high-resolution Ha and/or Paa images of the 
galaxies published by Alonso-Herrero et al. (2006 ) and adopted 
those as representative of the typical sizes of the actively star 
forming molecular gas complexes. 

In Fig.|6^ we have represented the derived KS law. The two 
surface densities, Esfr and Xdense, follow a coiTelation over more 
than 4 orders of magnitude in Sdense- An orthogonal regression 
fit to the data results in a KS-law of the dense molecular gas with 
a power index = 1 . 12 ± 0.04: 

logSsFR = (1.12 ± 0.04) logSdense +(-2.01 +0.12) (9) 

or SsFR- 0.010 S^^'l (10) 

The location of the new sample of LIRGs in the KS 
diagram of Fig. helps visualize that extreme starbursts 
(LIRGs/ULIRGs) and normal galaxies are not fully overlap- 
ping in this scatter plot. In order to quantify if a two-function 
power law significantly improves the overall fit we have split the 
sample into normal (Lir < 10^' Lq) and IR luminous galaxies 
(LiR > lO" Lq). FigurelSJ) shows the result of the two-function 
fit, which can be expressed as: 

logSsFR = (0.90 ±0.06) log Sdense +(-1.71 +0.15) (11) 

or2:sFR-0.02S0f;;^ (12) 
for galaxies with Lir < 10" Lq, and: 

logSsFR = (1.05 + 0.05)logSde„se + (-1.70 + 0.16) (13) 

orSsFR-0.02S^;fl (14) 

for local and high-z IR luminous galaxies with Lir > 10" Lq. 

Whereas the index of the power law is not substantially dif- 
ferent in normal galaxies (A^ = 0.90 + 0.06) and LIRGs/ULIRGs 
(A^ = 1.05 + 0.05) this fit reinforces the idea of duality hinted 
at by the solution that Gracia-Carpio et al. (2008) found using a 
smaller sample of LIRGs/ULIRGs. As illustrated in Fig.|6j3, the 
extrapolation of the KS law fitting normal galaxies to the high- 
est values of Sdense (~a fewlO'*M0pc"^) underpredicts the Zsfr 
in IR luminous galaxies by up to a factor ~5. This is a factor ~4 
larger than the total uncertainty on Esfr- Within the range of gas 
surface densities shared by normal galaxies and LIRGs/ULIRGs, 
Sdense~ 10^-2 IO-'Mqpc"^, the disagreement between the two 
laws stays within a factor 2 to 3, i.e., still a factor ~2-2.5 larger 
than the total uncertainty on Ssfr- 

We have evaluated the goodness of the two-function power 
law fit by a standard analysis. We find that decreases by 
a factor of 1 .4 in the dual fit compared to the single power law, 
an indication that the dual KS law for the dense gas qualifies as 
a better description of the data. This result is reminiscent of the 
bimodality found in the star formation laws derived in normal 
star forming galaxies and mergers (LIRGs, ULIRGs and SMGs) 
from CO line data (Genzel et al. lIOTOl Daddi et al. l20T0l l. We 
nevertheless analyse in Sect. |5]what are the potential biases of 
the two-function fit of Eqs. (11-14), motivated by the division 
of our sample at Lir = 10" Lq. 

In our derivation of KS laws we have used different methods 
to derive the typical sizes assumed to be common for the actively 
star forming regions and the dense molecular gas disks. This is 
certainly a source of uncertainty, akeady included in our eiTor 
budget. We nevertheless note that if there is any bias it will af- 
fect both the estimated Ssfr and Sdense to the same extent. These 



uncertainties would make all data points to be shifted along tra- 
jectories that coiTespond to straight lines of slope unity. The fits 
found for the KS laws are very close to the linear behavior and 
the described uncertainties associated to source sizes cannot thus 
create a fake segregation. 

5. Statistical assessment of the two-function fit 

The results found in the previous section suggest that the KS law 
for the dense gas in LIRGs/ULIRGs differs from that in normal 
star-forming galaxies. In this section we quantify this conclusion 
by addressing the statistical biases inherent to these results. 

5.1. Bimodality 

We first analyse if the star formation laws derived from HCN 
qualify as bimodal. The definition of bimodality is very strict 
as it requires the existence of two distinct populations within a 
given sample with little if any object connecting the two pop- 
ulations, reflecting two different modes. We can assess if the 
star formation laws derived from HCN qualify as strictly bi- 
modal by inspecting the distribution of depletion times. Fig. [T] 
shows histograms of fdep in the whole sample and in the groups 
of LIRGs/ULIRGs and of normal galaxies derived using dif- 
ferent approaches as to the assumed conversion factors. As 
expected from the difference in the mean values derived in 
Sect, m LIRGs/ULIRGs and normal galaxies tend to populate 
different regions of the joint histogram. This yields the impres- 
sion of a bimodal distribution with two unequally populated 
peaks. Assuming the non-universal conversion factors discussed 
in Sect. |6] certainly increases the spread in fdep between normal 
galaxies and LIRGs/ULIR Gs (Fi g.]?]). 

Hartigan & Hartigan ( 119851 ) proposed a test that measures 
the goodness of unimodality as the best description of a statisti- 
cal sample. The test measures the dip of a sample, defined as the 
maximum distance between a given empirical distribution and 
its corresponding best fitting unimodal distribution. For a given 
sample size, the deeper the dip, the less likely is that the distri- 
bution is unimodal. An updated version of the Hartigan dip test 
is implemented in the Rpackage, a free software environment 
for statistical computing Q. We have applied it to our sample of 
galaxies, in order to evaluate if the set of depletion times is com- 
patible with an unimodal distribution. The test indicates that the 
dip of an unimodal distribution could be as large as measured 
for the depletion times in ~9Q% of the cases, which otherwise 
implies that unimodality is far from being rejected. These admit- 
tedly high probabilities decrease to ~80% if we assume the re- 
vised conversion factors for Esfr and Edense discussed in Sect.|6] 
and could reach ~20% if we only changed the conversion factor 
for Sdense, 1-6., if wc adopted a similar approach to that of Genzel 
et al. dintob . 

We nevertheless conclude that we cannot claim that the star 
formation process is strictly bimodal on the basis of the deple- 
tion timescale distribution. This conclusion holds specially in the 
case where universal conversion factors are adopted. 

5.2. Universality of the KS law for the dense gas 

The fact that the two-function fit is a significant improvement 
over the one-function fit suggests nevertheless the existence of 
an underlying duality in the distribution. Due to our a priori ar- 
bitrary splitting of the galaxy sample into two groups, defined 

" http://www.r-project.org/ 
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Fig. 6. a) (left panel) Star formation rate surface density Esfr as a function of dense molecular gas surface density as traced by the 
HCN(l-O) line, Sdense, in different populations of galaxies for which we have an estimate of the molecular disk size. Symbols are 
as in Fig. |3] The solid black line is the orthogonal regression fit to the full sample of galaxies, b) (right panel) Same as a), but 
showing the two function power law fit to normal galaxies (grey line) and IR luminous galaxies (black line). Here, normal galaxies 
(LiR < 1O"L0) are represented by open squares and IR luminous galaxies (Lir > 1O"L0) by black filled squares as referred to in 
text. Errorbars show +30% uncertainties in Ssfr and 2dense- 




6.5 7.0 7.5 8.0 8.5 6.5 7.0 7.5 8.0 8.5 6.5 7.0 7.5 8.0 8.5 

log(tdop/yr) log(tdep/yr) logftdep/yr) 



Fig. 7. Histograms of depletion times for the dense gas in the whole sample of galaxies (solid line) and in the populations of 
LIRGs/ULIRGs (grey filled) and of normal galaxies (hatched). The depletion times in the left-hand panel are calculated assuming 
the same FlR-to-SFR and HCN-to-mass conversion factors for the two populations. In the middle panel, the two conversion factors 
are modified (see Sect.|6ll, whereas in the right-hand panel only the HCN-to-mass factor is changed. 



by the dividing line at Ljr = 10^' L0, the evaluation of average 
galaxy properties in the two sub-samples is subject to potential 
biases, however. The underlying reason is that the selection cut 
used to separate the populations is not orthogonal to the variables 
being fit. 

In the following we test the null hypothesis that a single KS 
law with a fixed power law index for both LIRGs/ULIRGs and 
normal galaxies (e.g., Eq. 9) is a suitable representation of our 
data. This hypothesis can be rejected if, on the assumption that 
it holds true, the measured values of three chosen observables 
are unlikely to result from pure random scatter. The observables 
used as constraints are: 



- The scatter of the data with respect to the one-function power 
law fitted to the data (Eq. 9): crigt = y/S S Rm/(N - 2). 
SS Rifit is the sum of squared residuals and N is the num- 
ber of sources. The observable crifn tells us the amount of 
random scatter we can plausibly assume. The value derived 
for (Tifit is 0.20 dex. 

- The difference, in median, between the logarithms of the 
depletion times of the two populations: Alog(tdep) = 
median{log(tdep))normai - median{log(tdep))(U)LiRG- The mea- 
sured value for Alog(tdep) is 0.45 dex. The parameter 
A log(tdep) measures the difference between LIRGs/ULIRGs 
and normal galaxies, although it is biased towards positive 
values. The reason is that Esfr is essentially proportional 
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Fig. 8. Results from testing the null hypotheses of a single KS law (left column) and of a double KS law (right column). Each panel 
shows the 2.5%, 25%, 50%, 75% and 97.5% percentiles of criflt (top), Alog(t(]ep) (middle) and S S Rxat/S S Rifn (bottom) as a function 
of o"rand- The measured values are indicated by dashed horizontal lines. A gray rectangle represents the 95% confidence interval for 
frand derived in each panel from the measured value. Dotted lines show the confidence intervals derived from (r\6t (top panels). 
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to LiR, so that a positive noise-driven error in the Lir of a 
galaxy simultaneously underestimates its tdep and increases 
the probability that the object is classified as LIRG/ULIRG. 
Superlinear KS laws as that of Eq. (9) (slope~ 1.12) enhance 
this bias. 

- The ratio of S S R2fit/S S Rm, where 55/?2fit is the sum of 
squared residuals when using two functions to fit the data as 
in Sect. 14.21 The measured value for this ratio is about 0.70. 
This ratio measures the improvement in the fit when separate 
KS laws are adopted for the two populations. 5'5'/?2fit/5'5'/?ifit 
is always smaller or equal than one, since a two function fit 
can always match or even improve the solution of the one 
function fit. 

In order to calculate the corresponding probability distribu- 
tions of the observables we proceed as follows: 

- We assume that the intrinsic distribution of data in the 
(logSdense,logZsFR) plane is given by the corrected values 
that result from the orthogonal least squares fit to our obser- 
vations derived in Sect. 14.21 These points follow Eq. (9) by 
construction. 

- We can evaluate the effect of the scatter of the distribution 
on the definition of the subsamples by scrambling our sample 
independently along the Sdense and Ssfr axes with a Gaussian 
scatter, hereafter denoted as tTiand. The galaxies are then re- 
classified either as LIRG/ULIRG or as normal after rescal- 
ing the IR luminosity of each object by the randomized-to- 
observed ratio of Esfr- This mimics the effect of errors on 
the classification. 

- The scatter parameter o-^and appears as a free parameter in our 
analysis, since there may be sources of random scatter other 
than the observational error (0. 13 dex). We take 10 values of 
o"rand from ~ dcx to 0.9 dex and perform 5x10^ randomiza- 
tions for each value. Each randomized set of data is analysed 
as the observed data to derive the probability distributions of 
our observables as a function of o-^^nd- 

The left column of Fig. [8] shows the 2.5%, 25%, 50%, 75% 
and 97.5% percentiles of criflt (top), Alog(tdep) (middle) and 
S S R2fiJS S Rifit (bottom) as a function of cTiand assuming the 
null hypothesis. From each panel we can derive a 95% confi- 
dence interval for cTiand (gray bands in the plots) by requiring 
that the probability to get the measured value lies within the 
2.5% - 97.5% range. As shown in the top panel, cr,and should 
lie in the range 0.18 - 0.24 dex to be compatible with our mea- 
surement of cTifit at the 95% confidence level. Outside that range, 
CTifit would be either smaller or larger than measured (0.20 dex) 
in 97.5% of the cases. Comparing the three panels, we clearly 
see that this confidence interval is formally incompatible with 
those derived from the measured Alog(tdep) (0.25 - 0.81 dex) 
andSSR2nt/SSRifit (> 0.42 dex). 

We can thus conclude that the universal one-function fit is 
not the best description of the data. Although the (Xsfr, ^dense) 
distribution does not qualify as strictly bimodal according to the 
dip test, the two-function fit to the data qualifies as a better de- 
scription of the star formation laws derived from HCN. None 
of the biases related to the intrinsic non-linearity of the law or 
to the a priori arbitrary cut in Lir to divide the sample are able 
to account for the observed differences between normal galaxies 
and LIRGs/ULIRGs discussed in Sects. l4TT]andl4;2] 

Furthermore, following a similar proof by contradiction ap- 
proach, we have carried out a similar test assuming a double 
KS law determined by the best two-function fit of Eqs. (1 1-14) 



discussed in Sect. 14.21 Equivalent plots to those discussed above 
are shown in the right column of Fig. [8] These plots clearly show 
that the hypothesis of a double KS law is formally in agreement 
with the measured values. Therefore, although the test cannot 
prove that the hypothesis of a double KS law is the only solu- 
tion, we can conclude that it is a better description of the star 
formation laws derived from HCN than a single KS law. We 
note that it could be possible to save the single KS law model 
by assuming that a significant fraction of the measured scatter is 
not random, but due to the existence of hidden variables relevant 
to the star-formation process like global dynamical time scales. 
This alternative is explored in Sect. 17.31 

6. Conversion factors and star formation laws 

We discuss in Sect. 16. li the validity of the Mdense - ^hcn(i-o) 
version factor adopted above for LIRGs/ULIRGs and consider in 
Sect. l6.2l the biases in estimating SFR for normal galaxies from 
^m^- consequences of revising conversion factors for the 
star formation laws obtained are discussed in Sect. 16.31 

6.1. The HCN-fo-Mdense conversion factor in LIRGs/ULIRGs 

It is currently accepted that the conversion factor between 
the CO luminosity (i'coii-O)-' mass of molecular gas 

(Mgas) is lower in mergers than in the Milky Way (MW) 
(q'^'^=4.8MqL'"'). The first conclusive evidence of a differ- 
ent conversion factor in mergers was found by Downes & 
Solomon ( I1998I I. These authors derived that the Mgas - ^co(i-O) 
conversion factor in ULIRGs is ~l/5 of the MW's value (this 
lowering factor is 1/10 in the case of Arp220). More generally, 
the analysis of CO line emission by radiative transfer models 
indicates that the physical conditions found in the CO clouds 
of mergers and nuclear starbursts (densities n(H2)-10^cm~^- 
lO'^cm"-' and brightness temperatures Tr^20-50 K) imply that 
the CO conversion factor, which scales as n(H2)'^^/TR, is typ- 
ically about (1/2 -1/4) X a^<^ (Solomon et al (19971 Scoville 
et al. I1997I I. A clear-cut argument in favor of a lower CO con- 
version factor in mergers is found in the fact that the canonical 
MW value of a'^'^ applied to ^^0(1-0) predicts gas masses that 
are larger than the dynamical masses in many LIRGs/ULIRGs at 
different redshift ranges (Tacconi et al. 120081 Daddi et al. 120101 
Genzel et aL l20T0l 

There have also been claims in the literature that the MW 
value of the Mdense - ^HCN(i-O) conversion factor, initially 
adopted in Sect. H) does not apply to mergers either (Gao & 
Solomon 2004a ; Gracia-Carpio et al. 2008). The high HCN to 
CO ratios measured in some ULIRGs/LIRGs (^0.1) make the 
requirement that Q'm?rgei-s-( V2-l/4)xa;^'^ mostly incompatible 
with the basic prescription that Mdense should not in all likeli- 
hood exceed O.SxMgas (Gao & Solomon l2004al l. 

The case of the ULIRG Arp 220 is paradigmatic in this 
respect. Downes & Solomon (1998) derived from the lower 
conversion factor for CO a total gas mass for the system 
Mgas~5xlO^MQ. The mass of dense gas derived from our HCN 
observations if we use the standard conversion version factor for 
HCN (aJJ^) is Mdense- l-2xlOi"M0, i.e., ~(2-3)xMgas. This re- 
sult indicates that a^™20 should be <(l/2-l/3)a^™ even in the 
extreme limit case where Mdense ^A^sas- The new data obtained 
in our work for other IR luminous galaxies corroborates this 
picture. Figure |9] shows the HCN to CO ratios as a function 
of Lir derived for the full sample of galaxies used in this pa- 
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per. This diagram is used to single out over-luminous HCN line 
targets in our sample, characterized by HCN/CO>0. 1 within the 
errors. About 50% of the LIRGs/ULIRGs in our sample are over- 
luminous HCN galaxies, identified by the red filled square sym- 
bols in Fig.|9l By contrast, only one normal galaxy in our sample 
is over-luminous in HCN. In these galaxies a^™ has to be low- 
ered by a similar factor as a^^ogrs relative to the MW canonical 
values. 



I I 1 1 1 1 1 I I I 1 1 III I I 1 1 1 1 1| I I I 1 1 1 ii| I ±1 I I 
1 — " - 




10" 1 



Fig. 9. The HCN to CO luminosity ratios as function of Lir for 
the full sample of galaxies analyzed in this work. Here, normal 
galaxies (Lir < 1O"L0) are represented by open squares and IR 
luminous galaxies (Lir > 1O"L0) by filled squares. Red mark- 
ers highlight those galaxies showing over-luminous HCN lines 
(HCN/CO>0.1). With the exception of one normal galaxy, all 
galaxies in this category are IR luminous. 



Following a similar approach to that employed to constrain 
the value of o'^a-fers' Gracia-Carpio et al. (2008 ) used the data 
obtained in a muftiUne survey of HCN and HCO+ of 17 LIRGs 
and ULIRGs to fit a^'^^ with Large Velocity Gradient (LVG) 
radiative transfer codes. The main outcome of Gracia-Carpio et 
al.'s study also backs up the claim that a^™ is ~3 times lower 
at high LpiR. Lower conversion factor for HCN in extreme star- 
bursts can be related to the particular chemical environement 
of the dense molecular gas in these sources where a hot-core 
like type of chemistry may prevail and enhance HCN abun- 
dances (e.g., Lintott et al. 2005). In addition, non-collisional ex- 
citation could also be responsible for the surprising strength of 
HCN lines in some IR luminous systems (e.g.; Aalto et al. 119951 
Garcfa-Burillo et aL l2006l Guelin et aL l2U07l Weil3 et al. l2007l l. 
In either case, the HCN conversion factor should be lowered. 

The independent lines of evidence mentioned above indicate 
that for a sizeable fraction of the IR luminous galaxies in our 
sample, the HCN conversion factor should be lowered. It re- 
mains to be proved that a similar correction should be applied 
to all the galaxies in our sample that qualify as IR luminous 
mergers. With this caveat in mind we re-derive the KS laws in 
Sect. 163] adopting a^'=^=l/3.2xa^ in LIRGs/ULIRGs. This 
is the same re-scaling factor globally adopted for mergers in the 
CO study of Genzel et al. (2010). We note that about 70% of 
the LIRGs/ULIRGs analyzed in our paper have Lir > 10"'^Lq, 
the limit beyond which mergers and strongly interacting systems 



start to dominate over less disturbed disk systems (Sanders & 
IshidaEOOll Alonso-HeiTero et al. l2009] l. 

6.2. The LpiR-to-SFR conversion factor in normal galaxies 

In Sect. |4]we used the IR luminosity of the galaxies as a proxy 
for their SFR. While this approach is valid for dusty galaxies, it 
has now become clear that in normal star forming galaxies one 
needs to account for both the obscured (traced by the IR lumi- 
nosity) and the unobscured SFRs (Perez-Gonzalez et al. 120061 
Calzetti et al. 2007; Kennicutt et al. 120091 ). In particular, com- 
binations of an IR tracer (e.g., the 24jt/m luminosity) with the 
observed H^^ luminosities are shown to work well. In this sec- 
tion we evaluate the corrections needed to account for all the 
SFR (obscured+unobscured) in our sample of galaxies. 

For the normal star-forming galaxies in our sample we 
made use of the SINGS sample data. We have eight galax- 
ies in common with the SINGS survey. We used the inte- 
grated observed (that is, not corrected for extinction) Ha and 
24jL/m luminosities given by Kennicutt et al. (2009) to com- 
pute the total SFR using the prescriptions given in their paper 
As described in Kennicutt et al. ( 120091 ). the use of total IR or 
monochromatic IR luminosities work both equally well to de- 
rive the total (obscured+unobscured) SFR. On average we find 
that SFR(24jum)/SFR(tot) = 0.5 - 0.6. This indicates that SFR 
based on an IR tracer underestimates the true SFR roughly by a 
factor of 2 in normal star forming galaxies (see also Kennicutt et 
al. |2009| and references therein). 

In LIRGs and ULIRGs, as mentioned above, the IR-based 
SFR indicators account for most of the total SFR. This was 
shown by recent integral field spectroscopy studies of a sample 
of more than 50 LIRGs and ULIRGs (Garci'a-Marm et aL l2059l 
Rodriguez-Zaurm et al. i20l iD . Therefore no corrections are ap- 
plied to derive SFR surface density for LIRGs and ULIRGs. 

We thus adopt in Sect. 16. 31 a correction factor of 2 to derive 
EsFR from in normal star forming galaxies. This is likely 
an upper limit to the true correction factor in normal galaxies, 
as we expect that the attenuation will vary more continuously as 
a function of the FIR luminosity. We note however that apply- 
ing a lower correction factor for normal galaxies would further 
reinforce the dual behavior of star formation laws discussed in 
Sect. I63](see Fig. 171). 

6.3. Revised star formation laws 

We have fitted a new two-function power law to the whole 
sample of galaxies using the revised values of Sdense (in 
LIRGs/ULIRGs) and 2sfr (in normal galaxies) derived accord- 
ing to Sects. 16.11 and 16.21 Figure [TOb visualizes the new star 
formation relation in normal galaxies and LIRGs/ULIRGs. The 
best orthogonal fit solution can be expressed as: 

logSsFR = (0.90 + 0.06) log Sdense + (-1-40 + 0.15) (15) 

orEsFR-0.04 2°f°^ (16) 

for galaxies with Lir < lO" Lq, and: 

logEsFR = (1.05 ± 0.05) log Edense + (-1-17 + 0.16) (17) 

orEsFR-0.07S^fl (18) 

for local and high-z IR luminous galaxies with Lir > lO" Lq. 

An inspection of Fig. [TOb indicates that, compared to 
Fig.|6]b, the dual behavior is significantly reinforced in the new 
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Fig. 10. a) {left panel) Same as Fig.|6]3j, but showing the revised two function power law fit to normal galaxies (grey line) and 
luminous infrared galaxies (black line) discussed in Sect. 16.31 Symbols and errorbars as in Fig. |6l b) {right panel) We show the 
revised depletion time scale as a function of L'^^q^ derived in normal and luminous infrared galaxies, as discussed in Sect. 16.31 
The dashed and continuous horizontal lines indicate, respectively, the average value of the depletion time scale in normal galaxies 
(tdep -50+5 Myrs) and LlRGs/ULlRGs (tdep ~14±1.4 Myrs). Symbols and eiTorbai-s as in Fig.|5] 



version of the KS laws. At the high end of Sdense values (~a 
few 10"* M0pc"^) the normal galaxy law underpredicts Zsfr in IR 
luminous galaxies by about an order of magnitude, a factor ~8 
larger than the typical uncertainty on Zsfr- Furthermore, within 
the range of gas surface densities shared by normal galaxies and 
LIRGs/ULIRGs, 2dense~2 10'-2 lO^Mopc ^, the factor 3 to 5 
disagreement between the two laws is statistically significant, 
being a factor ~3^ larger than the typical uncertainty on Ssfr- 
Instead, if we try to fit all the observations shown in Fig. [TOb 
with a one-parameter law over the whole range of IR luminosi- 
ties we obtain: 

logSsFR =(1.15 + 0.04)log2dense + (-1.60 + 0.06) (19) 

This solution has a 2.2 times higher than the two-function 
fit found above. 

Figure \Wb shows the revised depletion time scales of the 
dense molecular gas in normal galaxies and LIRGs/ULIRGs. A 
similar reinforcement of the duality in star formation laws can 
be identified if we compare the new estimates of tdep with those 
shown in Fig. |5] The value of tdep is on average a factor ~3- 
4 lower in LIRGs/ULIRGs (tdep -14+1.4 Myrs) compared to 
normal galaxies (tdep ~50+5 Myrs). This difference is a factor of 
2-3 larger than the total uncertainty of individual data points. 

As extensively discussed in Sect. |5j the revised star forma- 
tion laws of Eqs. (15-18), by increasing the spread between nor- 
mal galaxies and LIRGs/ULIRGs, make more plausible the use 
of the two-function fit as the best description of the data. 

7. Star formation models: local vs global 

The evidence of bimodality in star formation laws derived from 
CO line data is a challenge for star formation models relying 
only on local processes (e.g., Genzel et al. |2010l l. Our findings 
of a similar yet admittedly less extreme duality in the star for- 
mation laws derived from HCN, a better tracer of the dense 



molecular phase, and thus more directly related to star forma- 
tion, can put stringent constraints on models. By contrast with 
the highly scattered CO diagram, the smaller intrinsic disper- 
sion in the HCN star formation laws helps assess the verisimil- 
itude of duality with a higher statistical significance. We dis- 
cuss in Sect. 17.11 to what extent our results based on the HCN 
observations discussed in this paper can be used to test one of 
the key predictions of the star formation model of Krumholz 
&McKee (I2005I I: the constancy of the star formation rate per 
free-fall time (SFRff). We study in Sect. l7.2l the values predicted 
by Krumholz et al.'s models for the L^i^j L'^^^^y ^^^ luminosity 
ratios in different galaxy populations and use our data to bench- 
mark the models. Furthermore, we study in Sect. l7.3l if the inclu- 
sion of 'global' dynamical time scales offers a better description 
of the star formation laws derived for the dense gas. 

7.1. Observations vs models: is ttie star formatiori rate per 
free-fall time constant? 

Krumholz &McKee (120051) developed an analytic prescription 
for the star formation laws that are valid for different galactic 
environments. In their model star formation is globally very inef- 
ficient as it only takes place in a small subregion of a supersonic 
turbulent virialized molecular cloud. In this over-dense subre- 
gion, known as the cloud core(s), gravitational pull is able to 
take over kinetic energy and the collapse starts. If SFR is the 
total star formation rate in a galaxy, the star formation rate per 
free-fall time in objects of class X, SFRff _x, can be obtained as 
a function of SFR and the total mass of X objects, Mx, (or of 
their associated variables SFE or tdep), as well as of the fraction 
of star formation occurring in these objects, fx and the free-fall 
time, tff_x: 

fx SFR ?ff-x ?ff-x 
SFRff _x = ^ , " = fx SFEx %-x = /x (20) 
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Fig. 11. a) (/e/t panel) Star formation rate per free fall time (SFRff ; shown in the left Y axis) derived from the observed SFEjense in 
different populations of galaxies assuming a characteristic gas density for the HCN cloud nHCN(H2)=3xlO"icm-^. We show in the 
right Y axis the value of necNCHi) derived from SFEdense assuming a constant SFRff =0.02. This value lies within the most likely 
range for SFR)f~0. 011-0. 028 as determined by the star formation model of Krumholz & McKee (2005). We also highlight the 
average value of SFRff~0.0058 as determined by Krumholtz & Tan (I2007I I from a compilation of galactic and extragalactic HCN 
observations. Symbols are as in Fig. [3] Errorbars on SFRff (+42%) and Lfir (±30%) are shown, b) (right panel) Same as a) but 
obtained based on the revised values of Edense (in LIRGs/ULIRGs) and Esfr (in normal galaxies) discussed in Sect.|6] 



The star formation rate per free-fall time is a function of two 
parameters: the ratio of the kinetic energy (thermal plus turbu- 
lent) to the gravitational energy of the cloud (aviR) and the Mach 
number of the region (Ai) where Ai^ is roughly the ratio of the 
kinetic energy to the thermal energy. One of the key predictions 
of these models is that the dimensionless SFRff is expected to be 
low and lie within a quite restricted range <3% (~0.01 1-0.028) 
for the entire span of plausible values of ffyiR (~l-2) and Ai 
(~20-80). This offers a quantitative description for why star for- 
mation is so inefficient in a highly turbulent molecular medium. 

If we take HCN clouds as representative objects then 
SFEx=SFEdense=lAdep as defined in Sect. 14. II We can therefore 
estimate SFRff using Eq. |20] directly from observations, pro- 
vided that we assume a typical density for HCN clouds and 
suppose a value for fx- It is reasonable to assume that fx is 
about unity as all star formation in galaxies is likely taking place 
at densities exceeding the critical densities of the HCN(l-O) 
line (n^^(H2)~3xlO'*cm-Vr, where t is the opacity of the 1- 
line of HCN; e.g., Gao & Solomon l2004al) . The actual value 
of tff_HCN depends on the density of clouds emitting HCN. All 
in all, observations can be used to verify/falsify the prediction of 
Krumholz &McKee ( I2005l l's model about the behavior of SFRff. 

Krumholz & Tan ( 120071) used galactic and extragalactic ob- 
servations of a variety of objects to derive SFRff_x as a function 
of the characteristic ISM densities, which are probed by differ- 
ent molecular tracers (i.e., X-objects) in their compilation. A key 
outcome of their work is that SFRff _x seems to stay roughly con- 
stant with density. However, the results obtained from HCN sug- 
gest a significant deviation of SFRff_x from this general trend. 
The average value of SFRff derived from the HCN observations 
used in thek paper (SFRff_HCN ~0.0058) lies below the range 
predicted by Krumholz & McKee (I2005l )'s model. Krumholz & 
Tan ( I2007I I nevertheless assigned a large error bar to SFRff_HCN, 



implicitly assuming that the scattered distribution of SFRff de- 
rived from HCN data is mostly random^. 

We use below the HCN data discussed in this work to study 
the distribution of SFRff values derived in our sample of galax- 
ies. Our goal is to check if the new observations are compat- 
ible with the prediction of a constant SFRff. In this compar- 
ison we use the two versions of the star formation laws de- 
rived using the standard conversion factors {Case-I below) and 
the revised conversion factors discussed in Sect. |6](Cflie-// be- 
low). As discussed below, exploring different conversion factors 
is paramount if we are to probe the plausible range of values 
for SFRff, as these two quantities are interconnected through 
Eq. (20). 

7.1.1. Case-I: standard conversion factors 

If we assume that all HCN clouds have densities equal to 
n^" ^=3 xlO'*cm--^ (i.e., we take t=1) we can obtain SFRff from 
Eq. 1201 Figure [TTl ; represents SFRff as a function of Lfir ob- 
tained from the SFEdense values derived in Sect. H) i.e., prior to 
the correction of conversion factors discussed in Sect.|6] We note 
that the distribution of SFRff lies noticeably below the range pre- 
dicted by the model, as shown in Fig.fTTl;. In addition to this con- 
spicuous downward shift, the deduced SFRff distribution shows 
a systematic trend with LpfR and thus cannot be described as 
random. The order of magnitude increase in SFRff from nor- 
mal galaxies to LIRGs/ULIRGs is statistically significant and it 
is about a factor of 7 larger than the typical uncertainty of in- 
dividual data points (42%). The SFRff progression echoes the 
SFEdense trend with LpfR discussed in Sect. 14.11 



^ &umholz & Tan ( 12007b claim that uncertainties in observations 
translate into an order of magnitude errorbar. 
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To make compatible a constant SFRjf with the order of mag- 
nitude increase in SFEdense, the density of HCN clouds has to 
change notably from normal galaxies to LlRGs/ULlRGs, as 
nHCN(H2) scales as -SFE^^^^j,^. We represent in the right Y axis 
of Fig. \lTb the value of nHCN(H2) required to fit observations, 
assuming a fixed SFRff=0.02 for all galaxies, which is close 
to the average value predicted by the star formation model of 
Krumholz & McKee (2005). The required nHCN(H2) densities 
span two orders of magnitude from normal galaxies (^lO^cm^-*) 
to LIRGs/ULlRGs (~10'*cm"^). Of note, these densities are 
exceedingly low compared to nj^j^; this poses a problem for 
the standard scenario of collisional excitation of HCN lines. 
Extreme and thus unrealistic high opacities for HCN would be 
required to compensate for the low densities: thcn=3-300. As 
an additional relevant constraint, we would need to assume that 
opacities are a factor 10-100 higher in normal galaxies com- 
pared to LIRGs/ULlRGs. This requirement is at odds with the 
commonly measured higher molecular abundances of extreme 
starbursts ( Comb es [19911 Wild et al. [19921 Nguyen et al. [19921 
Krips et al. 120081 Gracia-Carpio et al. l2008l l. The current obser- 
vational evidence contradicts the existence of these abnormally 
low HCN densities, at least in the galaxies for which reliable es- 
timates of nHCN(H2) exist (e.g., Tacconi et al. |1994t Sternberg et 
al. [19941 Usero et al. [20041 Ki-ips et al. [20081 Gracia-Carpio et 
al.|2Ci08^. 

Based on the new HCN observations presented in this work 
and the use of standard conversion factors, we derive values 
for SFRff and/or nHCN(H2) that are well below the theoreti- 
cal expectations both for normal galaxies and LIRGs/ULlRGs. 
Furthermore the different densities required to fit the two popu- 
lations of galaxies with a common constant SFRff ~0.02 are far 
from the values derived in observations. 
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Fig. 12. We show the model predicted ratio of the star forma- 
tion rate or IR luminosity to the HCN(l-O) luminosity as a func- 
tion of the mean gas density n for the three cases examined by 
Krumholz & Thompson (120071): normal(dot-dashed line), inter- 
mediate(solid line) and starburst(dashed line) galaxies. The dot- 
dashed horizontal lines highlight the LFiR/iHCN(i-0) observed in 
normal galaxies and LIRGs/ULlRGs. The vertical dot lines in- 
dicate the location of the best agreement between observations 
and models as discussed in the text. 



7.1.2. Case-ll: revised conversion factors 

FiguredTJ? is similar to Fig.[TTl; but here SFRff values have been 
obtained from the revised estimates of the conversion factors dis- 
cussed in Sect. |6] for Sdense (in LIRGs/ULlRGs) and Ssfr (in 
normal galaxies). Not surprisingly, by adopting the revised con- 
version factors the global distribution of SFRff is shifted upward 
and, at the same time, the difference between normal galaxies 
and LIRGs/ULlRGs is increased. Most of the LIRGs/ULlRGs 
in Figure W\\ ) show now values of SFRj which lie within the 
range predicted by Krumholz & McKee (2005) with HCN den- 
sities very close to nj^^j^. The problem persists however in the fit 
of normal galaxies, where the predicted values for SFRff and/or 
nHCN(H2) are still well below the expected range. The paradigm 
of a common constant SFRff in all galaxies can only be saved if 
we allow this parameter to be as low as ~0.0035. In this scenario 
the HCN densities fitting the observed efficiencies would range 
from n^j^ in normal galaxies to ~10^cm"^ in LIRGs/ULlRGs. 
The fit in densities nevertheless requires to assume a value for 
SFRff that is a factor of 5-6 lower than expected by Krumholz & 
McKee' s models or derived from CO observations (Krumholz 
& Tan [20071) . 

We can thus conclude that the paradigm of a constant 
SFRff predicted by models within the range covered by the 
new HCN observations of normal galaxies and LIRGs/ULlRGs 
is not supported. The use of revised conversion factors allevi- 
ates the problem especially in LIRGs/ULlRGs. However, within 
the framework of local models it is not possible to fit the ob- 
served differences in the SFEdense between normal galaxies and 
LIRGs/ULlRGs using a common constant SFRff ~0.02 and a set 
of physically acceptable HCN densities. 
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Fig. 13. We represent log(tdep) versus log(tdyn) for normal galax- 
ies (LfR < 1O"L0, open symbols) and luminous infrared galax- 
ies (LiR > 10"Lo' filled symbols). 
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7.2. Observations vs models: the LFiR/incNci o) ^'Jf^if^osity 
ratio 

Based on the ISM model developed by Kmmholz & 
McKee (120051 1. Ki-umholz & Thompson (2007) predicted how 
the Lfir/L^^j^i^j luminosity ratios should change as a func- 
tion of the galaxy mean gas density h for three different object 
classes: normal galaxies, intermediate galaxies and starbursts. 
In Krumholz & Thompson's formulation, molecular clouds 
in normal galaxies are characterized by low kinetic tempera- 
tures (Tk-^Q K), as well as moderate-to-low Mach numbers 
(M=30), HCN abundances (X(HCN)=lxlO-'^) and HCN(l-O) 
opacities (thcn=0.5). These molecular cloud parameters are 
expected to be boosted to higher values in starburst galaxies: 
Tjf=50 K, M=80, X(H CN)^4 xlO'^ thcn=2 (See Table 1 of 
Krumholz & Thompson 120071 for details). The common enve- 
lope to normal and starburst systems is that molecular clouds are 
characterized by a roughly constant dimensionless star forma- 
tion efficiency measured by SFRjj . 

One key prediction of this model is that for mean densi- 
ties n<10"^cm"^, normal galaxies are expected to show higher 
SFEdense than galaxies classified as starbursts. Another key pre- 
diction is that the Lfir/L^j^j^^j luminosity ratio should remain 
roughly constant in a diagram that contains simultaneously data 
from normal galaxies (characterized by low n) and starburst 
galaxies (characterized by high h) up to densities «<10'*cm"-'. 
A different (less flat) appearance of SFEdense as a function of 
h can be foreseen at higher densities, however. Figure [121 a 
modified version of Fig. 2 of Krumholz & Thompson (2007), 
purposely explores the predicted SFEdense in the high den- 
sity regime: 10"^cm"^<n<10^cm"^. Following the convention of 
Krumholz & Thompson (120071 ) we represent the LFiR/iHCN(i-o) 
and SFR/L[jpj^|.j ratios as a function of the mean gas den- 
sity h for normal (dot-dashed fine), intermediate (solid line) and 
starburst (dashed fine) galaxiefl As argued by Krumholz & 
Thompson (2007"), the models for normal galaxies and starbursts 
are seen to converge in the high density regime («>10^cm"^). 

As discussed in Sect. 14.11 the observations presented 
in this paper show that the LFiR/iHCN(i-0) ^^^^'^ is not 
constant inside our combined sample: SFEdense is a fac- 
tor ~2-3 higher in LIRGs/ULIRGs compared to normal 
galaxies (<SFEdense>(LIRGs/ULIRGs) -1400+100 LqL''; 
<SFEdense>(normal)~6OO±7OL0L'"'). These data can be used 
to benchmark the predictions of the model shown in Figure [12] 
Two conclusions can be drawn from this comparison. First, in 
relative terms, the model is able to reproduce the increase in the 
SFEdense derived from observations as being due to a dramatic 
increase in the average gas densities. According to Fig. [12] 
if we assume that normal galaxies are characterized by h-2 
lO^cm^^*, extreme starbursts (LIRGs/ULIRGs) should be in a 
«=10^cm"^ regime to fit the observed change in SFEdense- In 
absolute terms, the model predicts values of SFEdense for normal 
galaxies and extreme starbursts that are marginally a factor ~ 1 .5 
above the observed valuefl This success of the model would 



This figure has been produced using the code of the 
model made publicly available by Mark Krumholz at 
|http : //w w w.ucolick. org/ ~ krumholz/ 

^ We note that the absolute scale along the right Y axis of the Fig. 2 of 
Krumholz & Thompson ( 20071 seems to be incorrect. A re-evaluation 
of this scale, derived from the output of the model, indicates that the 
predicted Lf\^l L'^^^^^ ^^ ratios are a factor of 2-3 lower than shown in 
the mentioned figure. Of particular note, this change helps finding a 
better agreement between observations and the model. 



be in apparent contradiction with the scenario described in the 
previous section, where we concluded that the star formation 
efficiency per free-fall time SFRff in the model, supposed to be 
~0.02, has to be lowered by a factor 5-6 to fit the observations 
(after the revision of conversion factors discussed in Sect. 17.1.21 ). 
As argued in the following there is no contradiction, however. 
Based on the output of the model of Krumholz & Thompson, we 
can actually estimate the implicit value of for the model 
molecular clouds at densities larger than nyjij^=3xl0'*cm"^. The 
HCN conversion factor depends, for each galaxy class, on h. 
Assuming the best fit values for h reported above, we derive that 
q;HCn_2_2^q^'-i jjj starbursts and -I-AMqL'^^ in normal 
galaxies in the model, i.e., a factor ~2-3 lower than the values 
of a^'^^ adopted in Sect. 16.11 Taking into account that the best 
compromise solution for LFiR/L[jf,j^^| implies ratios that are 
still ~50% too high (Fig.[T2l). we conclude that the model (with 
SFRff ~0. 02 and our preferred values for a^'--'^) is a factor of ~4 
off. This is similar to the factor of ~5-6 disagreement found in 
Sect. 17X21 

The conclusion common to the present section and Sect 17. II 
is that a rough agreement between the HCN observations pre- 
sented in this paper and the predictions of local models can only 
be reached if SFRif ~0.005-0.01 and/or if o'"™ is still a factor 
of ~a few lower than our favored values for normal galaxies and 
starbursts (LIRGs/ULIRGs). In qualitative agreement with this 
picture, the recent paper of Krumholz et al. (2011 ) has presented 
conclusive evidence that the best correspondence between the 
prediction of local models and the results of CO galaxy surveys 
requires lowering SFRff to 0.01, i.e., a factor of 1.5-3 below the 
value previously assumed by Krumholz & McKee (120051 ). 

7.3. Observations vs models: is 'global' dynamical time scale 
a key parameter? 

We have provided evidence that the conversion of dense molec- 
ular gas into stars is more efficient in LIRGs/ULIRGs than in 
normal galaxies. Although this has been ascribed to the exis- 
tence of very different local physical properties (densities) of 
the clouds in these two populations of galaxies, we have seen in 
Sects I7.ll and l7.2l that this hypothesis has difficulties in explain- 
ing the HCN results discussed in this paper (unless significantly 
lower SFRff and/or a^'^^ are assumed). Alternatively, the star 
formation rates could be regulated by the global large-scale dy- 
namics of the ISM, more perturbed in strongly interacting galax- 
ies (in our sample most of them are LIRGs/ULIRGs) than in 
normal galaxies. In this scenario we expect that the differences 
in SFEdense (or equivalently in tdep) between the two populations 
of galaxies can be reproduced by similar differences in tdyn, as 
SFEdense=l/tdep'^l/tdyn- To cxplorc this possibility, we compare 
in this section the distribution of dynamical time scales and de- 
pletion time scales of the dense molecular gas in our sample. 

Figure [T3] represents log(tdep) versus log(tdyn) for normal 
galaxies and LIRGs/ULIRGs. Depletion time scales of the dense 
gas tdep are derived as defined in Sect. 14. II Dynamical time scales 
tdyn are defined as a crossing time scale unit determined by the 
ratio between the radius of the star forming region, where the 
HCN emitting gas is expected to lie, and half the full width of the 
HCN line. When available, we have corrected tdyn for the value 
of the inclination of the galaxies. Whereas there is a clear segre- 
gation between the two populations of galaxies as to the values 
of tdep, the corresponding differences in tdyn are less significant. 

Figure [14] shows a histogram representation of the number 
density of galaxies as a function of tdep, tdyn and the tdep/tdyn ra- 
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Fig. 14. a) (upper panel) Histograms showing the number den- 
sity of galaxies (number of galaxies per unit interval, as shown) 
as a function of log(t(jep) for normal galaxies (Lir < 1O"L0, 
grey filled histogram) and luminous infrared galaxies (Lir > 
10"Lo, empty histogram), b) (middle panel) Same as a) but 
representing the number density of galaxies as a function of 
log(tdyn). c) (lower panel) Same as a) but representing the num- 
ber density of galaxies as a function of log(tdep/tdyn)- 

tio. We split our sample into normal galaxies (Lir < IO^Lq' 
grey filled histograms) and luminous infrared galaxies (Lir > 
1O"L0, empty histograms). Depletion times (top panel) roughly 
span an order of magnitude in each subsample. They peak 
around ~10 Myr in the case of LlRGs/ULlRGs and around 
~50 Myr in the case of normal galaxies. The clear separation 
between the two histograms makes apparent the higher star for- 
mation efficiency of LlRGs/ULIRGs. The middle panel shows 
the dynamical time scales of the studied galaxies. The two his- 
tograms extend across a common range of ~ 1 .5 orders of mag- 



log(ZsFR) = 0.8 X log(ldense) -0-5 X log(tdyn) + 2.6 




Fig. 15. Three-dimensional representation of the bi-parametric 
star formation law giving the star formation rate density Ssfr as 
a function of the dense molecular gas surface density Edense and 
the dynamical time tdyn. 

nitude, although the dynamical time scales in LIRGs/ULIRGs 
are about a factor of 2 lower on average. The existence of differ- 
ences between the two populations of galaxies can be expected, 
for many LIRGs/ULlRGs are mergers where the interaction has 
brought most of the molecular gas to smaller radii. We note how- 
ever that the differences reported for tdyn fall short of explaining 
the factor 5 difference in tdep. 

We represent in the lower panel of Fig. [14] the histograms of 
the tdepAdyn ratio. This is equivalent to the ratio between the mass 
of dense molecular gas and the mass of stars formed in one dy- 
namical time scale. Although the shape (intrinsic dispersion) of 
the two histograms is similar, tdepAdyn ratios in LIRGs/ULIRGs 
are typically a factor of 2 lower than in normal galaxies. Taken 
at face value, this suggests that the amount of dense gas con- 
verted into stars per unit of dynamical time is somewhat higher 
in LIRGs/ULIRGs. This result would contradict the idea that star 
formation is governed mainly by global dynamical time scales, 
as in this case we expect that the tdep/tdyn ratio is roughly con- 
stant, thus showing no significant differences between the two 
galaxy populations. 

The inability of 'global' star formation models to fully ex- 
plain the KS law derived from HCN, anticipated above, is bet- 
ter illustrated in Fig. [15] We show in this figure the three- 
dimensional representation of the bi-parametric star formation 
law that expresses the star formation rate density Xsfr as a func- 
tion of the dense molecular gas surface density Sdense and the 
dynamical time tdyn. The best fit to the HCN data is given by: 

logEsFR = (0.8+0.1) logEdense+(-0.5±0.2)logfdy„+(2.6+l)(21) 

While the power law index for Xdense (0.8±0.1) is not far 
from unity within the errors, the corresponding index for tdyn 
(-0.5±0.2) is not compatible with -1. This is against the theoret- 
ical expected power law index for £dense/tdyn ratio that should be 
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about ~1 (Kennicutt ll998l l. Furthermore, the goodness of the fit 
of Eq. |2T] estimated by a standard analysis, indicates a mea- 
ger 12% improvement with respect to the value ofx^ obtained in 
the single parameter fit of Eq.[T9] With the inclusion of 'global' 
dynamical time scales we are not able to significantly improve 
the fit of star formation laws derived from HCN. 

Several biases in the way the dynamical time scales are 
derived could explain why the tdepAdyn histograms shown in 
Fig. [14] are not coincident for the two populations of galaxies 
analyzed in this paper In particular, dynamical time scales in 
normal galaxies, where star formation is typically much more 
widely spread in the disks compared to LIRG/ULIRGs, could 
be underestimated. However LlRGs/ULIRGs do not behave as a 
monolithic family in this respect: some LIRGs, especially in the 
range Lir <10'^ '*Lo can have spread-out star formation (Hattori 
et al. 2004; Alonso-Herrero et al. 120061 Rodriguez -Zahurm et 
al. lSOTTI . 

8. Summary and conclusions 

We have used the IRAM 30m telescope to observe a sample of 
19 LIRGs in the 1-0 lines of CO, HCN and HCO+. The galaxies 
have been extracted from a sample of local LIRGs with available 
high-quality and high-resolution images obtained at optical, near 
and mid IR wavelengths, which probe the star formation activity. 
The new data presented in this work allow us to expand the num- 
ber of LIRGs studied in HCN-HCO^ lines by more than a factor 
3 when compared to previous works. The chosen LIRG sample 
has a range of HCN luminosities that overlaps to a significant 
extent with that of normal galaxies. This is a key to testing if 
LIRGs represent a transition point in star formation (SF) laws 
between normal galaxies and ULIRGs. In particular we study if 
the observed bimodality of star formation laws in galaxies de- 
rived from CO line data in previous surveys, can be extended to 
the higher density regime probed by HCN. We use these results 
to explore the validity of different theoretical prescriptions for 
star formation laws where local and/or global dynamical time 
scales are explicitly included. 

We summarize below the main results and conclusions of 
this work: 

- We find that the star formation efficiency of the 
dense gas (SFEdense) estimated by Lfir/L^j^j^^j lumi- 
nosity ratios is on average a factor ~2-3 higher in 
LIRGs/ULIRGs (-1400+100 LqL'"') compared to normal 
galaxies (-600+701© L'"'). With the addition of the new 
LIRG sample we confirm on a more solid statistical basis that 
the LpiR-i'HCNii-o) correlation is significantly super-linear 
(power law ~ 1.23+0. 05) when we include normal galaxies 
and LIRGs/ULIRGs in the fit. 

- We have derived KS laws relating the two surface densities, 
SsFR and Sdense in normal galaxies and LIRGs/ULIRGs. A 
fit to the full sample gives a power index - 1.12 + 0.04. 
However, LIRGs/ULIRGs and normal galaxies are not fully 
overlapping in this scatter plot. In order to quantify if a dual 
fit applies to the power law, we have split the sample into 
normal and IR luminous galaxies. While the (Xsfr, ^dense) 
distribution cannot be described as strictly bimodal, we find 
that a two-function power law fit with indexes close to unity 
qualifies as a much better description of the star formation re- 
lation for the dense gas. This result is seen to be solid against 
the statistical biases inherent to this analysis. 

- If we account for the different conversion factors for HCN 
in extreme starbursts and for the unobscured star formation 



rate in normal galaxies, for which we find new evidence in 
these observations, the duality in star formation laws is en- 
hanced. This result extends to the higher molecular densities 
probed by HCN lines the more extreme bimodal behavior 
of star formation laws, derived from CO molecular lines by 
two recent surveys (Daddi et al. 120101 Genzel et al. 1201 Oi l. 
The revised depletion time scales for the dense molecular 
gas show a significant difference between LIRGs/ULIRGs 
(tdep ~ 14+1. 4 Myrs) and normal galaxies (tdep ~50+5Myrs). 

- We have confronted the results of the new HCN observa- 
tions with the predictions of theoretical models in which the 
efficiency of star formation is determined by the ratio of a 
constant star formation rate per free-fall time (SFRjj) to the 
local free-fall time (tff) (Krumholz &McKee l2005l l. Within 
the framework of these models we find that it is possible 
to fit the observed differences in the SFEdense between nor- 
mal galaxies and LIRGs/ULIRGs using a common constant 
SFRjf and a set of physically acceptable HCN densities, but 
only if SFRff~0.005-0.01 and/or if a^'^^ is still a factor of 
~a few lower than our favored values. 

- We have studied if the inclusion of 'global' dynamical time 
scales offers a better description of the star formation rela- 
tions derived for the dense gas. This approach helps finding 
a marginally better universal relation for the dense molecular 
gas. In the best fit solution for EsFROcE|j^j^_,g t"^^, the power law 
index for Edense (n=0. 8+0.1) is not far from unity within the 
errors, i.e., in agreement with the theoretical expected value. 
Nevertheless, the corresponding index fortdy,, (m=-0.5+0.2) 
is not compatible with -1, a result whose theoretical impli- 
cations remain to be understood. 

A fit of the SF laws similar to the one discussed in Sect. 17.31 
was first explored by Genzel et al. (.2010) , who used CO line data 
in a sample of normal SF galaxies and extreme starbursts at dif- 
ferent redshift ranges. The two-parameter (Mg^s, tdy,,) SFR law 
derived by Genzel et al. (|20T0] | (SFRocM^^t t^^**) is successful at 
providing a universal SF relation. Of particular note, the power 
law indexes in their fitted SF relation are also found to devi- 
ate from the theoretical expected values (SFRocMj,^^ '^dyn-'' how- 
ever Altogether the SF relation obtained in Eq. [2T| from HCN 
data and the one obtained by Genzel et al. 62010) from CO are 
very similar to the KS laws obtained, respectively, by Gao & 
Solomon (2004b) and Kennicutt ( 1998) provided that we remain 
in the 5^sFR-2^(dense)gas projection of the SF laws. 

The bottom line, however, is that the different dynamical 
time scale ranges found for normal galaxies and extreme star- 
bursts are not sufficient to dissolve the duality observed in the 
SsFR-2(dense)gas projection of the SF laws (see also Krumholz et 
al. l201 11 1. A possible explanation for this failure can be searched 
for in the uncertainties that plague the definition and the obser- 
vational determination of tdyn. This is paramount if we are to 
compare systems that are known to be in very different dynam- 
ical states. So far the most widely accepted definition for tdyn, 
also used in this work, is equivalent to that of a crossing or orbit 
time scale for the system, which is determined from observa- 
tions. However, while normal galaxies can be described as dy- 
namically relaxed rotating disks, mergers are strongly interact- 
ing and dynamically very disturbed. In this context using a com- 
mon definition for tdyn is questionable. High-spatial resolution 
observations can help disentangle the role of large-scale dynam- 
ics in the small-scale spatially resolved SF laws both in normal 
galaxies and in extreme starbursts. 

Alternatively, the inclusion of global dynamical time scales 
might not be the key to estimating the correct scaling of the SF 
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relation and thus determining the absolute value of the SFR for 
different galaxy populations. From the observational point of 
view there is ample room for improvement in our understand- 
ing of how conversion factors for different tracers of molecular 
gas may differ in galaxies. As discussed above, this is a prereq- 
uisite to accurately calibrating SF relations and exploring their 
differences between galaxies. So far, the evidence for different 
conversion factors in normal galaxies and extreme starbursts is 
limited to CO and HCN lines. Multi-transition studies to be done 
for different molecular tracers will be paramount to tackling this 
issue. 

Besides providing a better handle of conversion factors, re- 
quired to convert intensities (I) into gas surface densities (Sgas), 
multi-line observations can yield an estimate of volume gas den- 
sities (pgas)- This will ease a more direct comparison between 
observations and models that predict how the SFR volume den- 
sity should depend on pgas, rather than on Sgas. In this con- 
text it is worthy of note that Krumholz et al. ( 1201 11 1 have re- 
cently argued that a simple, local, volumetric star formation law 
with SFRff ~0.01 can dissolve the apparent bimodality found by 
Genzel et al. (1201 Oi l if different scale-heights are assumed for 
the gas disks (hgas) in normal galaxies and extreme starbursts. 
While this is a plausible explanation, we note that the observa- 
tional constraints on how hgas might differ in different galaxy 
populations are currently scarce. The change in hgas required to 
make duality/bimodality dissapear in SF laws remains thus to be 
observationally validated. 
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